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Abstract 

(N 

We consider the low energy effective field model of graphene monolayer. Coulomb 
interactions are taken into account. The model is simulated numerically using the 
lattice discretization with staggered fermions. The two point fermionic Green functions 
are calculated. We find that in the insulator phase these Green functions almost do 
not depend on energy. This indicates that the effective field model (in its insulator 
phase) does not correspond to the real graphene. 



^ 1 Introduction 



It is well - known that without the Coulomb interactions the effective field model of graphene 
monolayer is a good approximation to the original tight - binding model. This effective field 
q ! model [U El El IU [HI El El E] operates with the continuum Dirac field living in the graphene 
sheet. This continuum model is used also when the Coulomb interactions are switched oiJE 
We suppose that it remains a good approximation to the tight - binding model when this 
effective field model remains in the semi - metal phase, i.e. it does not predict the appearance 
of the energy gap. However, as it will be explained below, we have some doubts that this 
model may be applicable for the small values of the substrate dielectric permittivity, where 
it predicts the appearance of the fermion condensate (see below). 

Recently, the effective field model of graphene monolayer with the Coulomb interactions 
taken into account was investigated numerically using nonperturbative lattice methods 0. 
The application of numerical lattice methods is justified by the fact that the Fermi velocity 
Vp is about 1/300. That's why the effective coupling constant a ~ 137^ i s large and the 
Coulomb interactions are strong. Therefore, nonperturbative effects may be strong. In 
H31 [IU [IHl EH HI ESI HI HSl EO] the effective low energy field model of graphene was 
investigated numerically using the lattice regularized model with staggered fermions 0. The 
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3 In this case the continuum Dirac field interacts with the dynamical field of the electric potential. 

4 Within the ranges of perturbation theory the effect of the Coulomb interaction on various physical 
quantities was investigated in a number of papers (see, for example, [51 UOj and references therein). 

5 Within the original tight - binding model the problem was considered analytically in ;T5] while in [21] it 
was investigated numerically. 



1 



main output of these investigations is that there exists the phase transition at a certain 
value of the effective coupling constant (3. This effective coupling constant is related to the 
dielectric permittivity e of substrate as follows |21j : 



(It is worth mentioning that due to the lattice artifacts this relation may be modified within 
a certain lattice realization of the effective field model of graphene. ) 

There is an evidence that this is the semi - metal - insulator phase transition. Namely, one 
of the possible condensates becomes nonzero at f3 < j5 c [HI [21]. In addition, the indications 
were found that the usual longitudinal conductivity vanishes at the position of the phase 
transition [21]. The possibility that the insulator phase may appear in graphene monolayer 
has also been discussed in another context (see, for example, [25J and references therein). 

The possibility that the effective low energy field model describes well the real graphene is 
not so obvious when the effective low energy model is in the insulator phase. Our conclusions 
are based on the direct measurement of the two - point Green function in the lattice regular- 
ized effective field model of graphene. The regularization is based on staggered fermions. We 
simulate the model using the same code that was used earlier by one of us during the work 
on the paper [21] . This code was tested in several ways (in particular, some previous results 
on the graphene monolayer [161 El EH 120] were reproduced). We demonstrate that in the 
insulator phase the Green function almost does not depend on energy while its dependence 
on the space - like momentum remains nontrivial. This means that the correlation time 
becomes negligible. At the same time the correlation length in physical units may remain 
nonzero or, even, infinite. Therefore, the physical energy of the fermion excitation tends to 
infinity and the given field - theoretical model is not self - consistent at the corresponding 
values of e and, thus, cannot describe the real physics. 

The paper is organized as follows. In Section [2] we consider the details of the model. In 
Section [3] numerical results are represented. In Section 0] we end with the conclusions. 

2 The effective field model of graphene monolayer in 
lattice regularization 

In the present paper we use notations adopted in [30] and [21]. The model contains two flavors 
(corresponding to spin) of the 4 - component spinors ip coupled to the electric potential A± 
(we work in the imaginary time representation). The Green function has to be considered 
in a certain gauge. The gauge freedom of the system corresponds to the transformation 

our numerical procedure we fix this gauge freedom via 
the condition A^x 4 , z) = for the 3D point z = 0. (We unfix the value of A 4 at a certain 
point on this line.) The two - point fermion Green function has the form: 



G = = \\ D^DA^ y e W (- 1 - J d 4 x[djA 4 ] 2 




a = 1,2; J, J 



1,2,3 



(2) 
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Figure 1: The values of go (left figure) and \g\ (right figure) at k\ — in the semi - metal 
phase ((3 = 0.2). The lattice size is 20 3 . Error bars are within 2% of the considered quantities. 
We denote here k t = A; 4 , k x = hi, k y = k 2 . 

When the model is considered in lattice regularization, the values of momenta belong to 
the Brillouin zone. The lattice regularization contains mass parameter m (for the details see 
[2T]). It has to remain nonzero for the numerical algorithm to stay at work. Physical results 
are to be obtained when the extrapolation to m = is made. 

In the lattice regularized model with staggered fermions the single Grassman variable \I/ 
is attached to the sites [31]. In terms of \I/ the free fermion action has the form: 

S = + £ £ [**a xi * x+t -* x+ .a xi * x ]), a xi = ( 3 ) 

x z i=l,...,4 

In order to return to the original spinor and flavor indices of the spinors the lattice is 
considered with even number of lattice spacings in each direction [31]. Let us subdivide this 
lattice into the blocks consisted of elementary cubes. Each block has 2 3 sites (two lattice sites 
in each direction). We denote the coordinates of the blocks by yi. Therefore, the coordinates 
of the lattice sites are Xi — 2yi + rji, rji — 0, 1. We define the new fields 

i%]: = lY,nrrnr a *2y +v (4) 

Here index a = 1, ...,4 is the spinor index while a = 1, ...,4 is the flavor index. Matrices $ 

have 4x4 components. But not all of these components are independent. Eq. (j3J) leads to 
the constraint 

rsi^r^ = % (5) 

This constraint reduces the number of flavors from 4 to 2. The free propagator of $ in 
momentum representation (of the blocked lattice) has the form (see [31] and also [30| I2T]): 



G 



(E F a ^-smk a - i(m - ^-(1- cosA; a )r5 <g> T 5 T a 
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Figure 2: The values of g (left figure) and \g\ (right figure) at ki — close to the phase 
transition (/3 = 0.08). The lattice size is 20 3 . Error bars are within 2% of the considered 
quantities. We denote k t = fc 4 , k x = k±, k y = k 2 . 

\ Ea r a sin k a + %{m - | £ a (l ~ cos k a )T 5 ® T 5 T a ) 

32[E a |(l-cosfc a ) + m 2 ] [ } 

Here Tj = Tj acts on the flavor indices while T matrices act on the Dirac indices. Momenta 
k are k x = ^ k 2 = S fc 4 = 2 ^ 4 2 + ^ , i^i, ^2, ^4 G Z; the lattice size is N 3 . At the 
end of the calculation one must set m = 0. The terms proportional to (1 — cos k a ) ~ k\ ~ a 2 
disappear in the continuum limit (the other terms are proportional to ~ a; here a is the 
lattice spacing). 

We suppose that when the interaction is switched on the form of the Green function is 
the same: 

G = g 1 T 1 + g 2 T 2 + g 3 T 4 + ig + ig a 5 T 5 ® T 5 T a . (7) 

The terms proportional to g£ are expected to be negligible in the continuum limit similar to 
the corresponding terms without the Coulomb interactions. The values of g a , a = 1, 2, 3 can 
be calculated as 

i z i y,z r\,r\' 

S(v'i - fa + 5 ia ]mod2)(G(2y + V ,2z + v ')) (8) 

Here (G(2y + rj, 2z + 77')) is the staggered fermion propagator in the external field averaged 
over the configurations of the U(l) gauge field and over the pseudofermion configurations. 
In addition we calculated the go component of the Green function as follows: 

9o(k) = WN 2 l N 2 N 2 E eik(Z ' V) E S(rf t - fa + <Mmod 2) (G(2y + V ,2z + 1/)) (9) 

1 2 t y,z rj,ri' 
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3 Numerical results 



We simulate the model at m = 0.01. We collected enough statistics to calculate the Green 
functions over all momentum space at j3 — 0.05,0.07,0.08,0.09,0.1,0.2 on the lattices 10 3 
and 20 3 . On the smaller lattice the Green functions are calculated using the direct inversion 
of matrices. On the larger lattice of size 20 3 we calculated the Green function using the 
stochastic estimators (for the description of the method see [21]). According to J2TJ [16] the 
values (3 = 0.05,0.07,0.08 belong to the insulator phase while the values (3 = 0.09,0.1,0.2 
belong to the semi - metal phase. We do not observe any qualitative dependence of the 
results on the lattice size. 

We analyze the data on the values of the Green functions and have found that there is the 
essential excess of \g\ = J 1 g\ + g\ + g\ at k a ~ (a = 1, 2, 3) over the average value within 
the momentum space lattice for = 0.2. The dependence of this quantity on momentum 
is represented in Fig. [TJ On this figure we represent the values of | <7 1 = \jg\ + g\ + g\ and 
go attached to the points of the dual lattice. The value at the point of the dual lattice is 
obtained via the averaging over the vertices of the corresponding cube of the original lattice. 
Namely, we plot the values \g\ = J ^ J2 a ,v 9a{v) 2 and g = | J2v9o{ v )i where the sum is over 
the vertices v of the given cube and over the components a = 1, 2, 3. On this figure the four 
peaks represent the single one due to the periodic boundary conditions. 

We observe that deep in the insulator phase (at (3 = 0.05) the Green function practically 
does not depend on fc 4 . Moreover, g 3 is negligible compared to g\ and g 2 ■ This means that 
different time spices correlate with each other very weakly, and that the system is described 
by the effective 2D model rather than by the 2 + 1 dimensional model. The dependence of the 
Green function on ki, k 2 demonstrates an essential excess of \g\ = \J 'gf + g\ at k\ = k 2 = 
over the rest of the momentum space lattice. The value of go at k\ = k 2 = in the insulator 
phase is essentially larger than the value of go in the semi - metal phase for fc 4 = k\ = k 2 = 
(see Fig. El again, the values are attached to the points of the dual lattice and are averaged 
over the vertices of the corresponding cubes). 

Close to the phase transition (j3 = 0.07,0.08,0.09,0.1) the maxima of \g\ and go as the 
functions of k^ at k\ = 0, k\ = 0, k 2 = are observed, in principle. However, the hights of 
these maxima are very small, while \g\ and go depend on k± very weakly (see Fig. [2]). That's 
why at m = 0.01 we observe smooth transition between the two regimes in the vicinity of 
the phase transition (its position is pointed out, for example, in [2U EE])- The first regime 
corresponds to the effective 2D description of the theory approached deep in the insulator 
phase. The second regime corresponds to the traditional semi - metal phase. 

4 Conclusions 

We have simulated the lattice regularized effective field model of graphene monolayer with 
the Coulomb interactions taken into account. We calculate the fermion Green function in 
momentum space. We consider several points on the phase diagram: deep in the the semi - 
metal phase, close to the position of the phase transition pointed out in [111 Ell [El E31 CEO 
[HI [20], and deep in the insulator phase. 
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At (3 = 0.05 (deep in the insulator phase) the Green function practically does not depend 
on /c 4 . This means that different time slices correlate very weakly. Moreover, the values of 
93(h) ~ (tfj(k)T4tfj(k)) are negligible compared with g%, $2 for any values of the momentum k 
(where g\{k) ~ {ip(k)Tiip(k)), gi ~ (ip(k)T2ip(k))). This means that deep in the insulator 
phase the energy of the fermionic excitation tends to infinity. Close to the position of 
the phase transition semi - metal - insulator we observe the intermediate behavior of the 
mentioned above quantities. Namely, there are very small hights of the peaks of | <7 1 = 

\J 9i + 92 + di an d #0 ~ (4>(k)4>(k)) as functions of fc 4 . 

This is confirmed also by the consideration of the results for the current - current corre- 
lator as a function of imaginary time represented in (21] (Eq. (16)). Namely, in Fig. 3 of 
[2T] the spectral density of this correlator is represented. It is clear from this Figure that in 
the insulator phase at e < 4 (i.e., for (3 < 0.8) the only maximum of the spectral density is 
at the frequencies /c 4 ~ 1/a, where a is the lattice spacing. For fc 4 << 1/a the values of the 
spectral density are much less than at fc 4 ~ 1/a. Therefore, the correlation time extracted 
from this correlator is of the order of the lattice spacing (in the insulator phase) . 

We consider the mentioned above results as an indication that at the sufficiently small 
values of the dielectric permittivity of the substrate the effective low energy field model does 
not correspond to the original tight - binding model. Therefore, it may have nothing to do 
with the reality. Most likely, here the discreteness of the graphene honeycomb lattice becomes 
important for the description of the physical phenomena and the excitations that are not 
described by the effective field theory play an important role. Therefore the conclusion of 
p31[]Il[]ia[Hl[22l[^ that there is the insulator phase of the graphene 

monolayer, seems to us questionable. 

It is worth mentioning that we measure our quantities at fixed m = 0.01 while the phase 
transition was observed in the behavior of the quantities extrapolated to m = 0. In order to 
make definite conclusions it is necessary to repeat the calculations described in the present 
paper for different values of m and to extrapolate the results to the value of m equal to zero. 
Also the dependence of the quantities on the lattice size has to be investigated. This should 
be a content of the further investigation. 

To conclude let us mention the recent work [23], where the experimantal results are 
presented with no sign of the insulator phase in graphene. 

The authors kindly acknowledge comments of M.I.Katsnelson, and a private communica- 
tion with G.E.Volovik, as well as the discussions with the members of the lattice ITEP group 
M.I.Polikarpov, P.V.Buividovich, V.I.Zakharov, O.V.Pavlovsky. The numerical simulations 
have been performed using the facilities of the supercomputer centers of Moscow University, 
Kurchatov Institute, and ITEP. This work was partly supported by RFBR grant 11-02-01227 
and by the Russian Ministry of Science and Education (program "Human Capital" and the 
contract No. 07.514.12.4028). 

6 This conclusion is made on the basis of the numerical investigation of the given effective theory. 
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Figure 3: The values of go (left figure) and \g\ (right figure) at k\ = in the insulator phase 
((3 = 0.05). The lattice size is 20 3 . Error bars are within 2% of the considered quantities. 
We denote here k t = k±, k x = hi, k y = k 2 . 
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